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ABSTRACT 

We use numerical simulations to investigate the importance of HD formation and cooling on the 
first generation of metal-free stars in a ACDM cosmology. We have implemented and tested non- 
equilibrium HD chemistry in an adaptive mesh refinement simulation code and applied it to two 
situations. (1) It is first applied to the formation of 10^ — 10® M© halos which form in the absence 
of any ionizing source ("unperturbed" halos). We show, in agreement with previous work, that HD 
cooling is of only marginal importance for most halos; however, we find that for the lowest mass halos, 
with masses a few times lO^M©, HD cooling can equal or surpass the H2 cooling rate. This leads to a 
population of stars formed in halos with effective HD cooling that are less massive by a factor of ~ 6 
compared to halos dominated by H2 cooling. (2) In the second part of the paper, we ionize the halos 
in order to explore the impact of HD cooling in the presence of an ample population of free electrons. 
This leads to cooler temperatures (due to the electron-catalyzed production of H2) implying somewhat 
lower resulting proto-stellar mass. Adding HD chemistry changes this by lowering the temperature 
further, to the level of the CMB. We find that HD cooling dominates over H2 cooling in the density 
range 10^ cm~^ to 10® cm~'^, but above this density, the temperature rises and H2 cooling dominates 
again. Because of this, the accretion rate on to the protostar is almost the same as in the H2 case (at 
least for accreted masses below 50-100 M©), therefore we argue that HD cooling in ionized halos will 
probably not result in a population of significantly lower mass stars. 

Subject headings: cosmology: theory - galaxies: high-redshift - galaxies:formation - stars:formation — 
methods: N-body simulations 
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1. INTRODUCTION 

The conditions surrounding the formation of the 
first stars out of primordial gas has been the focus 
of much work over the last few years, and ha s re- 
cently been summarized bv lBrornm" fc LarsonI ()2004D and 
ICiardi fc Ferraral ()2005f ). The picture that emerges is 
that of the first stars forming in halos with masses 
around 1 0® M© due to cool i ng fr o m molecular hydro- 
gen (e.g., [ Couchman fc Rees' '1986'; 'Haiman et al.! Il996l: 
Tegmark e t al. 1997; Abel et al. 2002; Yoshid a et all 



20031) . The actual masses of the first stars are less 



well determined, but in the last five years a number 
of groups have employed high-resolution hydrodynamic 
simulations, usmg both adaptive Euleri an codes (e.g., 
lAbel et al.ir2002l : IO'Shea fc Normanl|2007f ) and smoothed 
particle hy drodynamic s (SPH) L agrangian code s (e.g., 
Bromm et al. 2002; Yo shida et al. 2003. 2006; G ao et all 



2007() to answer this question. This work has generally 



resulted in the formation of a single, relatively high mass 
(of order IOOMq) star forming in each halo. The result- 
ing ionization front from the first star ionizes the origi- 
nal c l oud before any other stars can form (^V ha len et al.l 
[200l iKitavama et"an[200l I Abel et al.l[2007h . 

Much of the early work focused on single objects 
that formed within the most massive halo in the sim- 
ulation volume. More recently, some groups have ex- 
plored the parameter space by conducting multiple sim- 
ulations with different in itial conditions and box sizes 
(jO'Shea fc Normaiill2007f) an d with varying cosmological 
parameters ( Gao et al.l 120071 ). In general, a fairly con- 
sistent picture results: massive stars form in overdense 
regions (halo masses ~ lO^M©) less than a billion years 
after the big bang (10 z < 30). 



Most first star simulations have relied on a nine species 
chemical reaction network consisting of hydrogen, he- 
lium, molecular h ydrogen, and their associated ions (e.g., 
lAbel et all Il997| ). This simple model is possible be- 
cause big bang nucleosynthesis produced only lithium 
and lighter elements. While only trace amounts of H2 
are present initially within the primordial gas, it forms 
rapidly in the high density gas within collapsing halos, 
and once the density exceeds n ^ 10* cm"'^ the three- 
body formation process drives the fraction of molecular 
hydrogen to near unity. 

H2 has been shown to be effective in cooling primordial 
gas within dense halos, lowering the gas temperature to 
~ 200K. However, the H2 molecule itself is a relatively 
poor coolant since it lacks a permanent electric dipole 
moment. On the other hand, the HD molecule has a per- 
manent dipole, and cools more effectively at low temper- 
atures (< 200K). This has led several authors to suggest 
that HD will have an important effect on fi r st star forma- 
tion (|Puv et al.lll99l iGalh fc Pallal [l998l : iStancil et all 
[T998h . 

However, in low-mass primordial halos, the temper- 
ature never falls below about 200K, which is required 
for HD to dominate over H2 cooling, and so the gen- 
eral consensus has been that HD is unimp ortant for pri- 
mordi al, low-mass halos. For example, iBromm et al.l 
()2002|) included HD chemistry in 3D SPH simulations 
of first star formation and found that H2 cooling dom- 
inated : Yos hida et al. (2006) came to the same conclu- 
sion. iRipamontil (|2007l ) recently conducted a suite of 
ID simulations with full HD chemistry and found that 
for typical halos this is true, but that in low mass halos 
(< 3 X lO^M©) HD is formed in sufficient quantities to 
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be as important in cooling the gas as H2 . In one run, HD 
cooling lowered the temperature of the gas to ^ 707^, 3- 
4 times lower than when HD was neglected. iRipamontil 
|2007l) argued that the Jeans mass in such halos is low- 
ered by a factor of 10, implying that significantly lower 
mass stars will form in these regions. This is an intrigu- 
ing suggestion and one of the motivations of this paper is 
to test this idea using full three-dimensional simulations. 

While there is general (but not universal) agreement 
that HD is likely to be mostly unimportant in low-mass 
halos that have never been ionized, a number of au- 
thors have recently pointed out that HD will form in 
high abundance when the gas is ioniz ed (as originally 
pointed out bv iShapiro fc Kand |1987[ ) . For example, 
iNagakura fc Omukail (|2005l )" used one-dimensional mod- 
els to follow the evolution of HII regions inside massive 
halos. They found that for high enough masses, the gas 
would re-collapse and form significant amounts of HD 
due to the presence of free electrons. This cooling low- 
ered the gas temperature to the CMB temperature and 
so the authors sugge sted that low mass st ars (of order 
1 Mq) could result. iMachida et~all (|2005[ ) also looked 
at the impact of HD cooling i n the swept up shells of 
primordial supern ovae (see also ISalvaterra et al.ll2004f) . 

More generally, I Johnson &: Bromml ( 20061 ) used one 
zone models to show that HD cooling would be important 
whenever primordial gas was significantly ionized. This 
included halos forming behind strong shocks, in relic HH 
regions and in high-mass halos. They showed that the 
cooling was sufficient to quickly reach the CMB tempera- 
ture and argued that this w ould give rise to a Population 
H.51 (M ackev et all I2003D with masses 10 times lower 
than traditional Population HI stars. iGreif fc Bromml 
(|2006l ) went on to suggest that these stars generated the 
majority of the photons that ionized the universe. 

In this paper, we use three-dimensional simulations to 
explore the impact of HD cooling in two related areas 
of primordial star formation. The first is to investigate 
more carefully the role of HD cooling in "unperturbed" 
low-mass primordial star forming regions, paying par- 
ticularly close attention to the suggestion of Ripamontl 
(|2007l ) that HD may play a key role in the lowest mass 
halos. 

In the second part of this paper, we explore the impact 
that free electrons have on HD cooling in primordial ha- 
los. As noted above, ionization can occur for a variety 
of reasons, including relic HII regions (the kpc-sized re- 
gions ionized by nearby primordial stars) and shocks due 
to mergers with higher mass halos. 

We take a deliberately simple approach to the genera- 
tion of free electrons and we simply ionize and heat the 
gas instantaneously at a certain point in the evolution 
of our halos. This is an approximation to more sophis- 
ticated models which follow the radiative transfer of the 
ionizing radiation and its heating explicitly. Our method 

^ Stars which formed in primordial, metal-free gas that has been 
perturbed by a previous generation of stars have often been referred 
to as Population II. 5; i.e., an intermediate generation between Pop- 
ulation III and Population II. A recent proposal (O'Slicact al. 
120081) is to refer to the "first generation" stars formed in primordial 
gas as Population III.l, and the "second generation" stars formed 
in gas of primordial composition but affected by feedback effects 
from the first generation as Population III. 2, a nomenclature we 
will adopt in this paper. 



is very similar tolO'Shea et al.l (|2005D . As discussed more 
fully in lMesinger et al.l "( 2006f ). this approach has two dis- 
tinct differences compared to a more realistic treatment. 
First, the optically thin radiative cooling approximation 
means that the dense gas in the core may be improperly 
ionized (although if the source is in the center of halo 
itself, this is less of a concern). Second, the instanta- 
neous nature of the heating prevents the gas from being 
expelled from the halo. The advantage of our approach 
is that we are able to compare the effect of HD cooling 
with and without ionization in similar halos, and so we 
can see the key effects simply and clearly. Also, it is 
computationally faster and so we can study more cases. 
Finally, as we will argue later, many of our key results 
should still apply when carried out with more realistic 
treatment of the ionization (and indeed, do match other 
simulations where we can compare). 

In the next section, we describe the simulation code 
and our treatment and tests of HD chemistry. In H'S.ll 
we first examine the impact of HD on primordial halos 
without ionization, and then in ij3.2l we analyze the role 
of HD cooling in halos with ionization. Section |4] con- 
tains a discussion of the general effects of HD cooling in 
primordial star formation, including predictions for the 
final masses of stars based on the end state of the simula- 
tions. Finally, a brief summary of the results is presented 
in gSl 

2. METHODOLOGY 

To follow the formation and evolution of the primordial 
clouds from cosmological initial conditions, we use the 
Enzo adaptive mesh refinement (AMR) method. It is a 
well-tested cosmological hydrodynamics scheme that has 
previously been used to study a number o f astrophysical 
syste ms including clusters of galaxies (e.g.,lBrvan fc Void 
200l[ lYoungcr fc Bryan 2007), t he ISM CSlv z ct al] 
2005 ). the intergalactic medium (|Richter et al.l 120061 . 
and other applications. In particular, it has been widely 
used to study the forination of the first stars at high red- 
shift dA bel et al."2000', 2002; Machacck et al|l200li EOQI 
lAbel et al. 2007; O'Shca fc Norman 200 7.). 

The code is detail ed mo re fully in [ Bryan fc NormanI 
(fT99l . iBrvanI (fl999l) . and lO'Shea et al.l (I2004D . and"^ 



we only summarize it here. The code relies o n a hybrid 
particle- mesh solver (|Efstathiou et al.l Il985f) , using an 
N-body representation for the collisionless dark matter, 
and an Eulerian grid for the baryonic component. Enzo 
includes a sec ond-order accurate piecew i se parabolic 
meth od solver (jColella fc WoodwardI 11984 iBrvan et al.l 
11995) for solving the hydrodynamic equations. The res- 
olution of the grid is increased over the course of a simu- 
lation using the block -structured AMR technique (see 
,Berger fc Colellal Il989l) . On every level and at each 
timestep, the state of the grid is examined to determine 
if any regions have exceeded critical values in baryonic 
density, dark matter density, or Jeans length. If so, then 
the grid is refined in that region by adding more grid 
points. In this way, the most interesting regions of a col- 
lapsing medium receive the highest levels of refinement, 
and dynamic range large enough to correctly model the 
small-scale physics of cloud collapse in a cosmologically 
significant volume is achieved. 

As our method and application is v ery s imilar to 
lAbel et al.l (|2002D and iMachacek eFal] (|2003f) . we de- 
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Fig. 1. — Equilibrium abundances of various species obtained 
from steady-state runs of Enzo. Lines show the analytic expecta- 
tion from simple derivations based on the reaction rates (see Ap- 
pendix). Diamond symbols mark the values converged on by adi- 
abatic runs of Enzo(with all radiative cooling turned off), showing 
that the chemistry solver reproduces the expected values. 

scribe only the changes in detail^. In the next two sec- 
tions, we discuss first our enhanced chemical model and 
then the details of the simulation setup. 

2.1. Chemical Model 

The first stars formed in a chemically simplified envi- 
ronment consisting of lithium and lighter elements. Cool- 
ing from atomic hydrogen and helium is inefficient in 
the low temperature (< 2000K) gas within collapsing 
halos at z ^ 20, but should enough molec ular hydro- 
gen fo rm it can cool the gas to T ~ 200K. lAbel et ahl 
(|1997l ) and lAnninos et al.1 (|1997[ ) developed a chemical 
model for primordial gas following the evolution of the 
nine dominan t species H, H+, H ~, He, He"*", He++, H2, 
HJ, and e~. lAbel et aLl (|1997D demonstrated that the 
reaction rates for these species can be reduced to 19 col- 
lisional and 9 radiative processes and remain valid over 
a wide range of temperatures (1 < T < 10*K) and den- 
sities (n < lO'^cm"^). The rate equations are solved us- 
ing an implicit scheme based on the back ward differenc- 
ing for m ula (BDF) techniq ue described in lAnninos et al.l 
(|1997f) . lAbel et al.l (|2000fi modeled the formation of a 
primordial molecular cloud using this chemical model, 
following the cloud fragmentation and collapse until a 
density of 10''^cm^'^ and molecular hydrogen abundance 

^ The code used is nearly identi c al to the public version 
of Enzo described in lO'Shea et al.l l|2004l ) and available at 
I htt p : / / lea, ucsd .edu/software /enzo / [ 



of 10^'' was reached. In lAbel et all (|2002f) the forma- 
tion of molecular hydrogen via the three-body process 
was added, which extended the range of valid densities 
to n < 10^°cm~'^ and allowed the collapse to continue 
until the formation of a fully molecular, ^ IMq pro- 
tostar. No heating term for three-body II2 formation 
is included, an effect which can be important for the 
temperature evolution of the gas at n > 10^ cm~'^. In 
this paper we focus on an intermediate regime at densi- 
ties 10'^ — 10^ cm~'^ where the temperature is lowest. 
This regime is not sensitive to the high density core, and 
though we allow the simulations to continue to very high 
densities 10^^ cm~'^) we do not consider the core in 
our analysis. 

In this work, we introduce three additional species 
to the network: D, D"*", and HD. Reaction rates for 
the deuterium species are listed in Table [TJ A to- 
tal of eight reactions are considered for the deuterium 
sp ecies. Five of the six reactions of the minimal model 
of iCalh fc PaUa (T99E) ar e incl uded, using the updated 
rates from jGIalli & Pallj (|2002f ). We neglect photoion- 
ization of deuterium as we do not consider any exter- 
nal radiation fields other than the CMB, which does not 
contribute to this rate at z < 100. We include two ad - 
ditional deuterium reactions from 'Galh & Palla' (200?) 
that were not part of the Galli & Palla (1998) minimal 
model. Lastly, we include deuterium ionization by col- 
lisions with electrons, for which the rate is unknown 
and the equivalent hyd rogen rate is adopted (see, e.g., 
iClover fc ,T aDDsen"2007l). 

The cosmolog ic al deu terium abundance is taken from 
iBurles fc TviQ ((l998l) and is D/H = 3.4 x 10"^ (by 
nu mber). This is ~ 40% larger than the value obtained 
bv lRomano et all (|2003[ ) using WMAP first-year results, 
though it i s within t he range of other D/H measurements 
(sec Roma no et al.l[2003i . Table 1). 

We ado pt the r ecentl y computed HD cooling function 
of ILipovka et al.l ()2005[ ). Their rate calculation includes 
all radiative and collisional transitions for the J < 8 
rotational levels and the v = 0,1,2,3 vibrational lev- 
els. The rates are similar to previous formulations (e.g., 
iFlower et al.ll200"ol ) except at high densities and temper- 
atures (T > lOOOK) where the collisional ro- vibrational 
transi tions provide added cooling. From ILipovka et al.1 
(|2005l ). we use their equation 5 for the low-density limit, 
and in the LTE approximation we use the polynomial fit 
from their equation 4 with a density of nn = lO^cm"'^ 
(at which point the cooling function has saturated) . The 
high and low density limits are combined as 



A(nH,T) 



A(LTE) 



1 



A(LTE) 



(1) 



We prevent cooling below the CMB temperature by set- 
ting the HD cooling rate to zero when the temperature 
in a simulation cell falls below the CMB temperature at 

the current redshift 

FoU owing iBromm et al.l ()2002f ) and lJohnson fc Bromml 
(|200l . we verified the chemical network in Enzo by com- 
paring the equilibrium abundances of the species most 
important for cooling against values derived analytically. 
While Enzo is designed primarily for cosmological sim- 
ulations, it is possible to set up idealized problems and 
study the physics in detail. We created an adiabatic 
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Fig. 2. — Thermal evolution of an ionized gas with an initial 
temperature of 50000K under isobaric conditions with and without 
HD cooling (solid and dashed lines, respectively). Once a temper- 
ature of < 300K is reached, HD cooling dominates and the gas 
cools to essentially to zero temperature, whereas gas with H2 cool- 
ing does not cool below ~ 60K. Note that this is not a cosmo- 
logical simulation and heating from the CMB was not considered. 
The e volutionary phase s of the gas are very similar to Figure 2 
of Yos hida et al.l Il2007al) . though the transitions occur on slightly 
different timescales due to the different rates used. 

steady-state model consisting only of baryonic gas at 
an initial temperature and density uniformly distributed 
within a static (non-expanding) volume, and with radia- 
tive cooling disabled. Using this setup, we ran simula- 
tions ranging in temperature from lO'^K < T < lO^K 
and allowed each simulation to run until the abundances 
of all 12 species converged. The results are shown in 
Figure [11 the simulation data are in excellent agreement 
with the values derived analytically in the appendix. 

In addition to examining the equilibrium properties of 
an isothermal gas, we considered the temperature evolu- 
tion of hot gas that is cooling isobarically. This simplified 
model approximates the conditions of photoionized gas 
as may be expected for relic HII regions, SNe bubbles, 
or the shock ionized regions of colliding halos. Using 
the setup described above, but with radiative cooling en- 
abled and maintaining isobaric conditions by forcing nT 
to be a constant at each timestep, we followed the tem- 
perature evolution of the gas until a minimum tempera- 
tu re was reached . This test was motivated by the work 
of I Yoshida et~all (|2n07an : comparison of our Figure [2] to 
their Figure 2 shows that similar results are obtained by 
two different numerical models. Figure [5] shows that HD 
cooling is important in ionized gas, lowering the mini- 
mum temperature to lOK (no radiative heating from the 
CMB was considered in this run). 

2.2. Simulation Setup 

The framework for our simulations is a single box 1 
comoving Mpc on a side within a ACDM universe with 
cosmological parameters consistent with ISpergel et al.l 
(|200l : = 0.76, = 0.24, = 0.041, Hq = 
73 km/s/Mpc, and as = 0.74. The box is partitioned into 
a root grid of 128'^ elements and initialize d at z = 99 with 
128"^ dark matter particles using an .Eisenstein fc Hul 
(|1999f ) power spectrum with a spectral index of n = 0.95. 



The evolution of the dark matter particles is followed to 
2 = 10 using a low resolution run, then four of the most 
massive virial ized halos a re identified using the HOP al- 
gorithm (Eise nstein fc H ut 1998). 

Each of the four halos represents a potential site for 
primordial star formation. The particles identified with 
each halo at 2; = 10 are traced to their positions at 
z = 99, and three additional levels of nested static sub- 
grids are defined within the root grid such that all halo 
particles are contained within the highest resolution sub- 
grid at all times (the highest resolution region has an ef- 
fective resolution of 1024"^). The subgrids are centered on 
a given halo, so that while all the simulations share the 
same root grid and initial conditions, the high resolution 
subgrids are defined independently for each. The initial 
conditions are then regenerated to include both baryons 
and dark matter. In the highest resolution subgrid the 
dark matter particle mass is 27.2 Mq and the average 
baryonic mass in each cell is 5.6 Mq. 

In the high resolution runs this innermost subgrid is 
further refined up to a maximum of 25 levels, with each 
refinement increasing the cell resolution by a factor of 
two along each axis. Cells are flagged for refinement 
whenever the baryonic or dark matter mass in a cell is 
more than 4 times the resolution masses quoted above 
(Truelovc ct al. 1997). In addition, a cell is refined if 
the cell width Ax is less than 0.25 times the local Jeans 
length. The dark matter density is smoothed on a scale 
of 1 comoving parsec in order to avoid heating of the 
baryons due to disc reten ess effects of the dark matter 
particles, as in Abel et al. (2002:). When the simula- 
tions reach the maximum refinement level, the cell length 
within the highest refinement region is 46.1 AU and the 
gas mass resolution is 9.7 x 10^'^ M©. 

For each halo, two runs are performed, one including 
HD chemistry and one without, to isolate the effects of 
HD cooling and facilitate comparison with previous sim- 
ulations where HD was not included. We test the case of 
primordial star formation in ionized gas as discussed in 
the introduction by repeating the simulations for three 
of the halos while ionizing the full simulation volume at 
z = 20. 

3. RESULTS 

3.1. Unperturbed Primordial Gas 

The first case we consider is for unperturbed primordial 
gas (Population III.l), beginning at z = 99 and evolving 
normally until the formation of a protostar. For all eight 
runs (4 x H2 -I- 4 x HD) a core (proper) density of 
n ~ 10'^ cm~'^ is reached at z '--^ 15, indicating the col- 
lapse of a protostellar seed. The H2 mass fraction quickly 
climbs to 10"'^ at a radius of ~ 1 pc and effectively 
cools the gas, lowering the temperature to a minimum of 
~ 300K. As the collapse continues, the density exceeds 
n ~ 10* cm"'^ and the H2 cooling rate becomes indepen- 
dent of density. The temperature of the contracting gas 
rises, and once a density of 10^ cm~^ is reached H2 is 
formed rapidly through the three-body process until the 
core is fully molecular. The simulations are terminated 
when the core density exceeds n ^ 10^^ cm~'^, though 
we do not consider the inner core in detail and restrict 
our discussion to densities < 10^" cm""^ (see WI.l^ . The 
inner ~ 20 AU contains a fully molecular ~ IM0 proto- 
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Fig. 3. — Radial profiles of gas properties at the final timestep for simulations with unperturbed primordial gas. H2 runs are indicated 
by thin blue lines and HD runs by thick red lines. The two more massive halos, Halol and Halo2, show little difference between runs with 
and without deuterium chemistry. The two less massive halos, Halo3 and Halo4, show significantly lower temperatures within the collapse 
region, and are also less dense within the inner parsec. 



star, consistent with the resuhs of lAbel et alJ ()2002[ ). 
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Fig. 4. — Upper left: ratio of cooling time to crossing time for 
all 8 runs, with H2 runs in thin blue lines and HD runs in thick red 
lines. Halol is shown as solid lines, Halo2 as dotted lines, HaloS as 
dashed-dotted lines, and Halo4 as long dashed lines. The cooling 
time is much greater in the Halo3-HD and Halo4-HD runs where 
HD cooling is efficient due to the reduced cooling efficiency at low 
temperatures. Lower left: ratio of crossing time to dynamical time, 
which is very close to one for all cases and indicates dynamical 
stability. Upper right: mass fractions of H2 and H. The fraction 
of molecular hydrogen approaches unity within the innermost part 
of the cloud. Lower right: mass fraction of HD, which is enhanced 
by a factor of 20-50 in the Halo3-HD and Halo4-HD runs. 

Properties of the individual halos are given in Tabled 
The four halos span a factor of 2 in virial mass. At the 
virial scale, HD has little effect on the halo gas. The col- 



lapse to a protostar occurs at similar redshifts and virial 
temperatures whether HD is included or not. However, 
the properties of the halo gas well inside the virial radius, 
particularly in the region where the minimum tempera- 
ture is reached, are clearly affected for two of the four 
halos, as is evident in the radial temperature profiles of 
all four halos (Figure[3]). Within ~ 3 parsecs, the H2 and 
HD runs for HaloS and Halo4 diverge, with the HD runs 
having lower average densities and temperatures com- 
pared to the H2 runs. This indicates that HD cooling is 
effective in these halos, while for the other two halos it 
is not. 

Figure m shows that throughout most of the collapse re- 
gion, HD production closely follows H2 production with 
a factor of 10~^ reduction in number density. The excep- 
tion is a region of enhanced HD production in a spherical 
shell ~ Ipc from the core, corresponding to the mini- 
mum temperature. For Halol and Halo2, the HD/H2 
ratio reaches a maximum of - 2 X 10^4 in this region, 
and HD cooling contributes less than 20% of the over- 
all cooHng, which is governed by H2 (Figure E]). On the 
other hand, for Halo3 and Halo4 the maximum HD/H2 
ratio is ^ and HD provides more than 90% of the 
total cooling at densities of lO^cm"^ ^ n < lO^cm"'^. 
HD cools more effectively than H2 at low temperatures; 
when HD cooling is dominant the minimum gas temper- 
ature drops from ^ 300K to < lOOK. 

Interestingly, in the two HD-cooling halos the overall 
cooling is reduced, as can be seen in the ratio of the cool- 
ing time to crossing time (Figure 3]). Near the temper- 
ature minimum, the cooling time rises dramatically for 
the HD-cooling halos, as the gas temperature is near the 
CMB temperature and the cooling rate goes to zero. In- 
side this region, the temperature rises while the HD/H2 
ratio drops. Thus in the inner part of the cloud the gas 
remains cool (T > 250K) but there is less HD, so that 
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Fig. 5. — Contribution of HD to the total cooling in the HD 
runs. Line styles used for each halo are described Figure|4] The top 
panel shows temperature as a function of density. For reference, the 
results from H2 runs are shown with thin lines. The middle panel 
shows the contribution to the total cooling from HD, with thin lines 
showing the contribution of H2 cooling in the HD runs. The bottom 
panel shows the HD/H2 ratio, indicating where fractionation has 
occured. For Halol and Halo2, HD provides < 10% of the total 
cooling the gas, while for Halo3 and Halo4, the HD/H2 ratio is 
greater by 3 — 5 X and HD provides > 90% of the total cooling over 
a wide range in density. 

the overall cooling remains low (icooi/^cross ^ 20)'^. The 
decrease in gas temperature caused by HD cooling actu- 
ally reduces the total cooling rate, so that the gas will 
contract quasi-statically. 

An obvious difference is that the two halos with ef- 
fective HD cooling have the two lowest masses (some- 
what below 10^ Mq), while the two H2 coo ling halos 
have masses around 10^ Mq or more. iRipam onti (2007) 
noted that halos below ~ 3 x IO^Mq collapse more slowly 
because H2 is not produced in sufficient quantities to 
efficiently dissipate the heat generated by gravitational 
contraction. The slower collapse allows HD to build up 
until it becomes a substantial coolant . This is similar 
to the results of iBromm et al.l (|2002l ) , who simulated 
a 2 X IO^Mq halo with HD cooling that went through 
an extended phase of quasi-hydrostatic contraction and 
cooled to lower temperatures. In both of these studies, 
when more massive (> lO^M©) halos were tested with 
and without HD chemistry, HD was found to have lit- 
tle effect. We find that halos as massive as 8 x IO^Mq 
generate substantial HD cooling, but do n ot take longer 
to collapse. In both'Ri pamontil ()2007f ) and IBromm et all 
((2002;), idealized conditions of a single halo were consid- 
ered, whereas in our simulations the full merger history 
of a halo is treated self-consistently. Thus we must exam- 

^ The cooling time is defined as E/E and the crossing time as 
r/cs- Both are locally defined within a grid cell. 



ine the evolution of the gas in more detail to understand 
why HD cooling is effective in low mass halos. 

We find that HD cooling in unperturbed primordial 
gas is a highly sensitive process. Two effects combine to 
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Fig. 6. — Evolution of the HD/H2 abundance ratio relative 
to the cosmological D/H abundance ratio during the collapse (the 
maximum enhancement within the collapse region is plotted). The 
X-axis represents the core density, and each halo is followed once 
a core density of 1 cm~'^ is reached. Linestyles are as in Figure |4] 
Filled symbols correspond to individual output times, with the 
color representing the minimum temperature in the collapse region 
at that time. In equilibrium, the HD enhancement would be ~ 6 at 
250K and ~ 100 at lOOK. The strong HD enhancement in two of the 
halos only occurs after the core density has exceeded ~ 10* cm~^ 
and the temperature has remained < 200K. For HaloS, the HD 
enhancement increases from 6.5 at point A to 15 at point B, a time 
span of 0.7 Myr. Halo4 undergoes an increase in HD enhancement 
from 5.5 (point C) to 14 (point D) within 0.3 Myr. While HaloS 
continues to increase in HE) enhancement, for both halos this brief 
period of strong enhancement signals the point at which HD cooling 
goes from marginal to significant importance. 

drive the effectiveness of HD cooling in low temperature 
primordial gas. First, the primary formation channel for 
HD, 

D+ + H2 H+ + HD, (2) 

is endothermic by 462K, so that HD formation is favored 
over H2 by a factor e^^^/T (jGlovedl2007h . Second, the 
HD cooling rate per molecule is greater than that of H2 
below a critical temperature. This critical temperature 
is determined by the gas density and the abundances 
of H2 and HD, and is generally i n the range 150-250K 
(see the discussion and Figure 1 in lGloveil[2007h . Should 
H2 cool the gas to the critical temperature, HD cooling 
will become important and can cool the gas even further. 
This results in a runaway effect, as HD is formed rapidly 
in the cool gas. 

Indeed, we do not see a slow buildup of HD as in 
iRipamontil (|2007t l: rather, the HD enhancement occurs 
rapidly in the latter stages of the collapse. Figure [H] 
shows the HD enhancement in each of the four halos, 
expressed in terms of the HD/H2 abundance ratio di- 
vided by the cosmological D/H abundance ratio (both 
values are ratios of number densities, the cosmological 
D/H value is stated in i i2.ip . At an early stage in the 
collapse, when the core density is ^ lcm~'^, the maxi- 
mum HD enhancement is ~ 8 for all halos. When HaloS 
reaches point A, the minimum temperature is 220K, and 
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Fig. 7. — Distribution of the gas in the n — T plane for two halos, 
at a time when the central density is ~ 5 X 10^ cm"''. Contours 
represent fractions of (0.1, 2.5, 10, 100) X 10~^ of the halo virial gas 
mass, which is ~ W^Mq for Halo2 and ~ 7x IO^Mq for Halo3. The 
upper line shows the critical temperature at which the HD cooling 
rate per unit volume is equal to the H2 cooling rate, calculated 
at the corresponding density and temperature and assuming that 
riHj = 10"'' nil S'lid njjD = 3 X lO"** njjj , which are typical values 
for the gas in these regions (see Figure[5ll. The lower line shows the 
the temperature at which HD provides 5 times the cooling of H2. 
The upper panels show the gas distribution for both the H2 (left 
panel) and HD (right panel) runs of Halo2. Most of the gas is above 
the critical region; this halo does not exhibit effective HD cooling. 
The lower panels show that for HaloS the cool gas is shifted to 
higher densities and a much greater amount of gas is within the 
critical region; this halo does exhibit effective HD cooling. 

the HD enhancement is still only 6.5. Within 0.7 Myr, 
point B is reached, where the minimum temperature has 
dropped to 150K, and the HD enhancement has jumped 
to 15. Halo4 undergoes a similar process, at point C the 
minimum temperature is 225K and the HD enhancement 
is 5.5, but 0.3 Myr later at point D the temperature has 
dropped to lOOK and the HD enhancement is 14. Both 
halos collapse from ricorc = lcm~'^ to ricoro — lO^'^cm^'^ 
in ~ 100 Myr, thus most of the HD enhancement occurs 
in less than 1% of the total collapse time. In fact, we 
find that HD cooling exhibits a 'critical' nature, in the 
sense that either HD cooling is very important to the 
cooling or it is not important at all. This behavior stems 
from the rapid HD fractionation and the runaway effect 
of cooling. 

The conditions for effective HD cooling are then that 
the gas must cool by H2 cooling to the critical tempera- 
ture and remain there long enough for a buildup of HD 
to occur. This is most likely to occur in lower mass ha- 
los. Halos with large gravitational potentials have higher 
virial temperatures and more dynamical activity. Merg- 
ers from infalling gas clumps will disrupt pockets of cool 
gas. Lower mass halos, which collapse more uniformly, 
are better able to form dense pockets of cool gas where 
HD can form. 

We also find evidence for the critical nature of HD 
cooling by looking in detail at the low temperature gas 
in our simulations. Figure [7| shows the n — T plane for 
Halo2 and Halo3 at a time when the central density is 
5 X 10^ cm~'^, just before the large HD enhancement seen 
in two of the halos (Figure [6]). It is clear that the lower 
mass halo, HaloS, has much more gas at higher densities 
10''^cm"^) and low temperatures (T < 200K) than 



Halo2. This gas is well below the critical temperature 
and is cooling effectively by HD. While some gas in the 
more massive Halo2 is just below the critical line, most 
of it remains above and HD cooling is not effective for 
this halo. Our r esul ts, combined with t he findings of 
iRipamontil (|2007t l and lBromm et al.l (|2002t ). suggest that 
a set of conditions must be met for effective HD cooling 
in unperturbed primordial gas, and these conditions are 
most likely to occur in < lO^M© halos. 

3.2. Primordial Star Formation in Ionized Regions 

The second case we consider is primordial star forma- 
tion in ionized regions (Population III. 2). In sites of pri- 
mordial star formation, molecular hydrogen forms pri- 
marily through the H^ channel, 

H + e-^H-+7 (3) 
H + H- ^ H2 + c- (4) 

a reaction which is catalyzed by free electrons and thus 
enhanced in ionized gas. Ionized primordial gas could 
be found in shocks from colliding massive halos, or in 
the photoionized regions surrounding the first genera- 
tion of stars. The latter case is of particular interest, 
as the end state of > lOOM© primordial stars is prob- 
ably either direct collapse to black holes or complete 
destr uction in pair instability supernovae (jHeger et al.l 
I2OO3D . In either case a sizable HII region will result, 
and as the first stars themselves form in large overden- 
sities these fossil HII regions will host dense gas clouds 
that could be sites for subsequent star formation. As 
H2 forms rapidly in the ionized regions, the gas cools 
to T^ lOOK, where HD formation is highly favored and 
HD cooling can drive the temperature even lower, essen- 
tially to the CMB temperature. Th e resulting low Jeans 
masses imply lower mass sta rs (Naga kura fc Omukail 
l2005HJohnson fc Bromn]|2006l ). 

We explore the formation of primordial stars in ionized 
gas through the simple approach of ionizing the entire 
simulation volume aX z — 20. This technique has been 
similarly utilized with Enzo by to'Shea et al. (2005) but 
without HD chemistry. While this method is not equiv- 
alent to a full treatment of radiative transfer, it does 
allow us to examine in detail the effect of HD cooling in 
ionized regions for otherwise identical halos, while being 
much less expensive computationally. 

We used three of the simulation volumes described in 
§ 12.21 so that the evolution to z = 20 is identical to the 
unperturbed case (§ 13. ip . but the late stage evolution 
begins with a high abundance of free electrons. As pre- 
viously, the simulations are performed once without HD 
chemistry and once with HD included. 

Table|3]shows the properties of the halos for the ionized 
runs. Compared to the unperturbed case (Tabled, the 
collapse for a given halo occurs somewhat earlier and at 
lower virial masses and temperatures in ionized regions. 
The resulting density, temperature and enclosed mass 
profiles are shown in Figure [H As in Figure [H the HD 
simulations show distinctly lower temperatures than the 
H2 runs in the range 0.01 - 10 pc, but now the effect is 
independent of mass. 

In Figure [9] it is clear that the molecular hydrogen 
fraction is greatly enhanced over the unperturbed case 
as expected. The extra H2 lowers the gas temperature 
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Fig. 8. — The same as Figure |3] but for the flash-ionized runs. 
Similar results are found for the three halos examined, in all cases 
temperature minimum approaches the CMB temperature at the 
redshift of collapse. 

to ~ 130K even without HD cooling, a factor of > 2 
lower than achieved in the H2 cooling unperturbed ha- 
los. When HD chemistry is included, it forms rapidly 
in the low temperature gas. This can be seen in the 




10" 10"' 10" 
radius [pc] 

Fig. 9. — The same as Figure |4] but for the flash-ionized runs. 
The molecular hydrogen abundance is > 10~^ throughout due to 
the enhancement of the H~ channel by the free electrons. The 
cloud is dynamically stable (tcross/idyn ~ !)• HD is greatly en- 
hanced within Ipc, where T ~ lOOK. 

large bump in the HD fraction at ^ Ipc, which reaches 
a maximum of 4 x 10^^, nearly two orders of magnitude 
greater than seen in the unperturbed halos. As a result, 
the temperature is driven to the CMB floor, and HD pro- 
vides > 95% of the total cooling in the low temperature 
zone (Figure [TU)) . 

It is interesting to compare the radial development of 
the H2 and HD abundances shown in Figure [H Mov- 
ing inward from a radius lOOpc, the H2 mass fraction 
rises to a maximum of ^ 6 x 10"'^, and remains essen- 



tially constant until a radius of ^ lO^'^pc. On the other 
hand, the HD mass fraction steadily rises from ^ lOOpc 
to Ipc, as the temperature drops from ~ lOOOK to 
~ lOOK. Over the next decade in radius, the tempera- 
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Fig. 10. — The same as Figure [5l but for the flash-ionized runs. 
Note that HD cooling is now important for all halos, accounting 
for 95-99.7% of the cooling in the T < lOOK region. 

ture is roughly constant near the minimum, and the HD 
mass fraction remains at its peak of ~ 3 x 10~^. In 
this region HD fractionation is nearly complete, and the 
HD/H2 ratio reaches 4 x lO"'^ (Figure fTO)) . At a density 
of ^ lO^cm^'^, the HD cooling function saturates and 
the temperature rises, leading to a sharp drop in the HD 
abundance, while the H2 abundance is still roughly con- 
stant. In the inner O.Olpc, H2 forms rapidly through the 
three-body process while the temperature again climbs 
to > TOOK, so that the HD abundance again follows the 
H2 abundance reduced by a factor of 7 x 10~^. 

4. DISCUSSION 
4.1. General effects of HD cooling 

In general, HD cooling is only important once H2 cool- 
ing has lowered the gas temperature to ~ 200K. This only 
occurs once the gas is already quite dense (n ~ lO'^cm^'^). 
Thus HD cooling does not affect which halos will form 
stars, but it does alter the properties of the gas in the 
region of the temperature minimum. 

For halos with effective HD cooling, the minimum tem- 
perature is not only lowered by a factor of a few com- 
pared to runs without HD, the location of the minimum 
is shifted to a smaller radius and hence higher densities. 
When H2 cooling dominates the minimum temperature 
occurs at a density of n lO'^ — lO^cm"'^. When HD 
cooling dominates, the minimum temperature occurs at 
a density of n '--^ 10* — lO^cm"-^. 
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This effect is easily understood by considering the 
properties of the H2 and HD coohng functions. First, 
it is useful to assume that the gas is thermally balanced 
with heating due to gravitational collapse balancing ra- 
diative cooling. The gravitational energy input rate per 
unit volume can be approximated as 



form in less massive halo s, cool by HP, and have masses 
< IOOMq. According to iHeger et"an (|2003f ). metal-free 



fGM(r) 

n 

V r 



(5) 



To derive this simplified expression, we have assumed 
that the enclosed mass M(r) is proportional to r, or al- 
ternately that p{r) oc r^^, as can be seen in Figures [3] 
and [51 The dynamical time idyn oc n"^/^. 

Over a short range in temperature, the cooling function 
can be roughly approximated as 



A oc 



(n < ricr) 



(6) 



depending on whether the density n is above or below 
the critical density Ucr- For both the H2 and HD cool- 
ing curves, a is quite large, typically of order 3-5 in the 
interesting range. By equating these two expressions for 
the heating and cooling we can derive a relation between 
density and temperature. When we do this, we see that 
the temperature reaches a minimum at the critical den- 
sity, where it transitions from T j^-O-S/q ^^qX ^ n"-^/". 
The critical density for H2 is n ^ 1 0'^ — lO^cm^'^, while 
for HD it is n ~ 10^ - 10^ cm'^ (|Lipovka et all 120051 : 
iFlower et al.ll20dol ). This explains the shift in the loca- 
tion of minimum temperature to higher densities as seen 
in Figures [H] and [TO] (and Tables [5| and [3]). In addition, 
the high value of a ~ 4 explains the relatively flat tem- 
perature profiles seen in the simulations. The increase in 
the temperature at densities beyond the critical density 
ultimately leads to a shutoff of HD cooling and, as we 
see in the next section, important implications for the 
masses of the resulting stars. 

4.2. Masses of the First Stars 

We have shown that HD cooling can be as important 
as H2 cooling for some low mass primordial gas clouds 
and in ionized regions. We can now examine whether 
HD cooling changes the nature of the collapse enough to 
alter the mass function of the first stars. Figure [TTl shows 
the accretion timescales for all simul ations, ca lculated 
according to the Shu isothermal model (|Shulll977[ ). where 
m = O.975C5/G' (we obtain similar but noisier results 
using the instantaneous mass infiow rate, m = Anr'^Vrp). 

For primordial halos, the accretion times are similar 
for all of the H2 cooling halos, suggesting ~ 3OOM0 
stars will form after 10^ years of accretion. However, 
the two HD cooling halos accrete significantly less mass 
for a given accretion time. This results from the de- 
creases in temperature and enclosed mass at r ~ Ipc 
when HD cooling is effective. This is an important scale 
for star formation as this is the region at which the 
Bonner-Ebert mass is at a minimum. Our results sug- 
gest that stars of ~ 5OM0 will form, much lower than 
is typically seen in simulations of first generatio n star 
formation. We conclude, in agreement with Ripa montil 
(|2007( ). that there are two mass scales for the first stars. 
Some (most) stars will form in halos of IO^Mq, cool by 
H2, and have masses > lOOM©. The remainder will 



-p 

c 
o 

-H 

-P 

U 
O 
O 



unperturbed 




1000 



Fig. 11. — Accretion times onto protostars in all simulations 
calculated using the Shu rate, with t^cc = m/rh. Linestyles are 
the same as Figure |4] The upper panel shows results from the 
unperturbed runs, while the lower panel shows results from the 
ionized runs. In the unperturbed case, HD has little effect on the 
expected final mass for the stars in Halol and Halo2, while having 
a strong effect for Halo3 and Halo4. The accreted mass is reduced 
by a factor of ~ 6. In the ionized case, HD cooling reduces the 
accreted mass by < 2. The mass scale for primordial halos with 
effective HD cooling is quite similar to that of ionized halos. 

stars with masses > IGGM© will produce pair-instability 
supernovae, while those with masses 40 — IOOMq will col- 
lapse directly to black holes. Thus HD cooling may be 
important in setting the mass scale of the first stars and 
in the later evolution of the interstellar medium. 

The inner regions of the ionized halos are quite sim- 
ilar in density and temperature, but when HD cooling 
is included the temperature is lower by a factor of > 2. 
The accreted mass is affected by a similar amount; the 
lower panel of Figure [11] shows that HD cooling lowers 
the expected mass of primordial stars in ionized regions 
from > IOOMq to ~ 60Mq. This mass scale agrees well 
with the results of lYoshida et al.l (|2007lJ ). who exam- 
ined the development of a protostar in a primordial HII 
region with SPH simulations that included HD chem- 
istry. They followed the evolution of the protostar to a 
central density of lO^^cm"'^ including a detailed proto- 
stellar evolution model that traced the evolution to the 
zero-age main sequence. They predicted that a ~ 4OM0 
star would form in the HH region, more massive than 
earlier simulations had predicted. Our finding of a simi- 
lar mass scale for three different halos suggests that while 
second generation primordial stars are not as massive as 
the first generation, they are still quite massive. 

We note that extrapolating a final stellar mass from 
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estimations of the accretion rate is highly uncertain and 
some caution should be taken in intrepreting the stellar 
masses. It is unclear whether a direct relationship exists 
between accretion rate and final stellar mass, whether the 
accretion is spherical (Omukai & Palla 2003) or through 
an a ccretion disk (Tan fc McK ee 20M; McKee fc Ta.nl 
l2007f) . Taking these uncertainties into consideration, it is 
clear that HD cooling has a significant effect on both the 
temperature and density evolution of lower mass Popu- 
lation III.l halos. The total mass within a given radius 
is much lower when HD cooling is effective (Figure [3]), 
implying stars of much lower mass will form. In ionized 
halos, the differences are much smaller (Figure [S]), im- 
plying that in this case the resulting stars will not be 
significantly less massive. 

We do not see signs of fragmentation in projections 
of the inner collapse regions. We also do not see evi- 
dence for the formation of rotationally-supported disks, 
as I'circ/i'kcp < 1 in all our simulations. This is not a 
surprise, as we do not expect HD cooling to cause frag- 
mentation ~ like H2 it has a steep cooling function and so 
is stable to thermal instabilities. In addition, as we see in 
Figure [21 the cooling time never drops below the dynam- 
ical time, another important criterion for fragmentation. 

4.3. Rate uncertainties 

The results presented here depend on the accuracy of 
the chemical model. We showed in § 12.11 that the chem- 
istry solver numerically reproduces the expected abun- 
dances in equilibrium; however, the calculat ed values are 
only a s acc urate as the r ates given as input. iGlover et al.l 
(|2006f) and lGloveil (|2007f) have noted that two of the rates 
most strongly ti ed to H2 form ation are also the most un- 
certain (see also IGlover fc A bel 2008). The first, the H" 
channel, is less relevant to the unperturbed case but may 
be important in ionized regions. This uncertainty will 
propagate into simulations with HD, as HD pr i maril y 
forms through H2. As discussed in IGlover et al.l (|2006f ). 
this can be an important source of uncertainty in these 
calculations, and we reiterate the need for improved rate 
measurements and calculations. 

The second uncertainty discussed bv lGloveil ()2007l ) is in 
the three-body H2 formation process. This uncertainty 
will not qualitatively affect our results, as the three-body 
process is only important in the inner core of the col- 
lapse region, where the temperature climbs up to > 500K 
and HD is less important. We did not include an anal- 
ogous three-body process for HD formation, though for 
the same reason we do not expect this to be an important 
process. 

Another source of uncertainty is in the cosmological 
D/H ratio. As discussed in § 12.11 there is a factor ^ 2 
uncertainty in cosmological measurements of this value. 
Figure [6] shows that HD cooling halos undergo a rapid 
fractionation phase; this is unlikely to be sensitive to 
small changes in the cosmological deuterium abundance. 
In addition, the final amount of HD enhancement differs 
by a factor of 2 for the HD cooling halos, yet the fraction 
of cooling from HD is similar for both cases (Figure [6]) . 

5. CONCLUSIONS 

In this paper, we have implemented a minimal chemical 
network for tracking the formation of HD in the Enzo cos- 
mological hydrodynamics code. We have verified the net- 



work using both equilibrium and non-equilibrium tests 
and then applied this to the formation of primordial 
stars in high-redshift halos with cosmological initial con- 
ditions. We examined a range of halos under two differ- 
ent assumptions for the initial electron fraction. In the 
first case, the halos were unperturbed, with an initial 
electron fraction as predicted by standard cosmological 
evolution in the absence of ionizing backgrounds. As in 
previous work, the first collapsed objects form in 10^^^ 
M0 halos at 2: ~ 15. In the second case, we ionized the 
halos at 2; = 20 to investigate the impact of free electrons 
on the collapse. Although this simple 'flash' ionization is 
not a realistic description of, say, a nearby photoionizing 
source, it does allow us to investigate the impact of HD 
cooling in the presence of a source of free electrons. Our 
primary results are as follows. 

• For Population III.l star formation in high mass 
halos (roughly 10^ Mq and above), HD coohng is 
unimportant. However, for the lower mass halos in 
our sample there was a factor of 10-100 increase in 
the HD fraction (compared to the high mass ha- 
los). This, is turn, caused HD cooling to dominate 
H2 cooling over the density range 10 cm~'^ to 10'' 
cm~'^ (corresponding approximately to 10~^ pc to 
1 pc, in radii). The HD cooling permitted gas to 
reach nearly down to the CMB temperature. 

• We found that in order for HD cooling to be im- 
portant, a sufficient amount of gas needed to be 
at a low enough temperature and a sufficiently 
high density for HD to form preferentially over 
H2. However, once this critical line was crossed, 
rapid HD fractionation coupled with runaway cool- 
ing quickly led HD to dominate the cooling rate. 
This requirement for lower temperatures helps to 
explain why low mass halos formed HD preferen- 
tially over high mass halos. 

• For the ionized case (Population HI. 2), HD cooling 
was important in all halos examined regardless of 
mass. HD cooling was efficient enough for gas to 
reach the CMB temperature over an even larger 
range of radii. 

• In all cases, once the density increased above the 
HD critical density (around lO'^"^ cm""^), the tem- 
perature began to rise again and invariably in- 
creased by about ~ 500 K, leading to the suppres- 
sion of HD relative to H2. This means that the 
interior solution (r < 10~^ pc) relies only on H2 
cooling and so looks quite similar to the no-HD 
case. In addition, no sign of fragmentation was 
found. 

• Despite the fact that HD cooling led to lower tem- 
peratures (often down to the CMB temperature), 
we found that the impact on the predicted mass 
accretion rates was relatively mild. It was largest 
for the unperturbed case, where the low mass halos 
were predicted to form stars a factor of six times 
lower than the no-HD case (from BOOMq to ^ 50 
Mq, assuming an accretion time of 10^ years - 
see SI for a discussion of the uncertanties in ac- 
cretion mass estimates). For the ionized case, the 
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predicted accretion rates were nearly identical with 
and without HD, leading to at most a factor of two 
reduction in the mass of the predicted stars form- 
ing out of ionized regions (from ~ 100 Mq to ~ 50 
Mo). 

These results have a number of interesting implica- 
tions. First, they suggest a way for at least some rel- 
atively lower mass stars to form in primordial (or near 
primordial) halos, as some interpretations of the ob- 
served extremely l ow metallicity halo s tars have required 
(jTumlinsonI l2007t iKomiva et all l2007l ) . However, they 
also indicate that HD cooling by itself will not form sub- 
stantially lower mass stars (e.g. by a factor of ten) as 
some previous analytic work has suggested. In particu- 
lar, we have not found any sign of a new class of Pop- 
ulation III. 2 objects (at least for the initial conditions 
surveyed). 

Our work suggest two avenues of future research. The 
first is a more careful and systematic investigation of the 
conditions for HD cooling in "unperturbed" halos. While 
our results indicate that HD cooling is more important 
in low mass halos, it is not clear how this depends on 
merger histories or stochastic effects. More work is re- 
quired to find the fraction of halos which form via the HD 
channel for a complete sample of halos. In the ionized 
case, our methodology for introducing free electrons was 
deliberately simple. It would be use ful to examine more 
realistic ionization m odalities (e.g.. IWise &: Abell2007t 
lYoshida et aIll2007aD . 
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APPENDIX 

In this appendix we show the derivation of the analytical formulas used to test the equilibrium abundances of the 
chemical network, as described in § I2.1 | and shown in Figure [1] The chemical network for all reactions not involving 
deuterium is taken from the work of I Abel et al.l ()1997f ). The complete set of reactions and their associated rates are 
given in that work, for brevity we focus here on the species most relevant to cooling in the primordial gas. For ease 
of comparison, rates listed in this appendix are also denoted R#, with the number corresponding to the rate number 
used in Abel ct al. (1997) (and also within the Enzo source code). 

The equilibrium abundance of ionized hydrogen is calculated from the coUisional balance of H+ (|Abel et al.l ()1997f ) 
rates R1-R2): 

H + e-^H+ + 2e-, (1) 
R+ +e- ^R + hv. (2) 

By assuming ng = 71^+ , we derive 

n-ti+ = p-nu- (3) 
The abundance of neutral hydrogen is found by assuming nu tot — + '^h+ • 

The two processes for creating molecular hydrogen include the dominant channel and charge exchange with 
(rates R8 and RIO), 

H + H-^H2 + e^, (4) 
H++H^H2+H+. (5) 

The three processes governing the destruction of H2 are the reverse of the above charge exchange reaction, as well as 
collisional dissociations by electrons and neutral hydrogen (R11-R13), 

H2+H+^H++H, (6) 
H2 + e-^2H + e^, (7) 
H2 + H^3H. (8) 

Thus the equilibrium abundance of H2 is given by 

hnnnii- + kionB_n„+ 

= I n — TT^- (9) 

fcll»T-H+ + ki2ne + fcl3?T-H 

Although the H2 abundance in the low-temperature gas (T < lO'^K) typical of primordial star formation is dominated 
by the H~ channel, we do not neglect in our equilibrium calculation as it is important at higher temperatures. In 
addition to reaction ^ (also denoted Rll), the following two reactions produce (R9 and R17): 



H + H+ ^ H+ + hi^, (10) 

H"+H+^H+ + e". (11) 

is destroyed by reaction (O (also denoted RIO) and the following reactions (R18, R19): 

H+ + e"^2H, (12) 

H++H-^H2+H. (13) 

The equilibrium abundances for H2 and are found by solving equation ^ simultaneously with 

fc9nHnH+ + fcll?^H+'^H2 + fci7"H-?^H+ /, .X 

n„+ = ; . (14) 

H~ is an important intermediary in the formation of H2 . It is formed primarily by photo-attachment of H and 
(R7), 

VL + e- ^R- +hu. (15) 

H~ is destroyed in a series of collisional processes, including reaction (fTTj) (also denoted R17) and the following rates 
(R14-R16): 

H-+H^H + 2e-, (16) 

H"+H^2H + e~, (17) 

H"+H+^2H. (18) 
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This leads to the equihbrium equation 

n^- = ^^^""^ . (19) 

fesTT-H + kiiUe + ki^Tin + kienY[+ + fci7nH+ 

From Table [1] is primarily created by coUisional ionization and destroyed by recombination and collision with 
neutral hydrogen, 

"D+ = 1 TT ■ (20) 

HD is created and destroyed in four exchange reactions involving H2 (Table [1]), 

k52nu+nu2 + fc54'^D'^H2 ,01 ^ 

'^HD = T TT • (21) 

We note that the rate for coUisional destruction of H2 by H (rate ([8]), or R13) used in Enzo is significantly higher than 
the rate for the same process used bv lBromm et alJ ()2002l ) for the temperature range 10"^ < T < 10^. The equil i briuin 
abundance of H2 at T ~ 15000 K in our calculation is lower by a factor of 6 than that of I Johnson fc Bromml (|2006l 
see their Figure 1). Simulations of primordial star formation do not reach such high temperatures so this difference is 
unimportant to the work presented here. 
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TABLE 1 

Reaction rates for deuterium species 







Reaction 




Rate Coefficient 
(cm'^ 


Reference 


Dl 


D + 


e — » 


D+ ^ 


- 2e- 


See expression in reference 


A97/1'' 


D2 


D+ 


+ e- - 


D 4 


- 7 


See expression in reference 


A97/2'' 


D50 


D + 


H+ ^ 


D+ - 


f H 


2.00 X io-iOTO-402e-37.i/T _ 3 31 ^ lo-"Ti-48 


GP02/5 


D51 


D+ 


+ H 


D + 


H+ 


2.06 X 10-10rO-396e-33.0/T + 2.03 X 10-9T-0-332 


GP02/6'' 


D52 


D+ 


+ Ha - 


H+ 


+ HD 


10-9 X [0.417 + 0.846 log T-0.137(logT)2] 


GP02/2 


D53 


HD 


+ H+ - 


^ Ha 


+ D+ 


1.1 X lo-Se-^s/T 


GP02/4 


D54 


D + 


Ha ^ 


H + 


HD 


1.69 X 10-10e(-''680/T+198800/T2) 


GP02/1<= 


D55 


HD 


+ H ^ 


■ Ha H 


- D 


5.25 X io-iie(-*«o/T+i73900/T2) 


GP02/3'^ 



Note. — The references are I Abel et"aI|[T997l (A97) and lGalli fc Palla >'200? (GP02), with the rate 
numbers used in those papers shown after the reference. The first column gives the rate numbers 
according to the system used internally within the Enzo code. 

The coUisional ionization and recombination rates for deuterium are taken from the equivalent 
rates for hydroge n.^ Rate 6 fr om [G alli & Palla 2002 contains a typo and was obtained from the 
original reference. ISavin l l2002l .'^ IRipa monti (2007) noted that these rates become unphysically large 
at T < lOOK, following the example of that paper we set the rates at T < lOOK to be the same as 
the rate at T = lOOK. 



TABLE 2 

Properties of the simulation halos in the unperturbed runs 



Run 






Rvir 




T 

rmn 








(IO^Mq) 


(pc) 


(K) 


(K) 


(lO^cm-S) 


Halol-H2 


14.55 


10.84 


213.0 


1610 


268 


5.5 


Halol-HD 


14.71 


10.24 


207.0 


1566 


286 


4.1 


Halo2-H2 


15.35 


9.62 


194.7 


1563 


322 


3.6 


Halo2-HD 


15.34 


9.82 


196.2 


1585 


300 


3.2 


Halo3-H2 


14.88 


6.72 


177.9 


1195 


290 


31.6 


Halo3-HD 


14.74 


6.94 


181.5 


1211 


66 


23.7 


Halo4-H2 


14.56 


9.37 


202.9 


1462 


268 


5.0 


Halo4-HD 


15.02 


7.93 


186.4 


1347 


103 


85.2 



Note. — The runs are labeled "H2" to indicate a run without 
deuterium chemistry and "HD" when deuterium is included; they 
are otherwise identical. The redshift of collapse is determined as the 
redshift when the central density exceeds lO^'^cm" the virial mass, 
radius, and temperature are calculated at this time as well. The 
minimum temperature Tmin is calculated at the final timestep and 
found by fitting a cubic polynomial to the n-T plane in logarithmic 
space, in the region around the temperature minimum (e.g., the 
upper panel of Figure [Sjl. The value of Tmin and corresponding 
density n{Tmin) are determined from the minimum of the fitting 
function, and illustrate the effect of HD cooling in the region of 
lowest temperature. 



TABLE 3 

Properties of the simulation halos in the ionized runs 



Run 










T ■ 








(IO^Mq) 


(pc) 


(K) 


(K) 


(lO^cm-^) 


Halol-H2 


13.53 


12.72 


240.5 


1673 


132 


1.3 


Halol-HD 


15.06 


7.38 


181.5 


1286 


55 


9.1 


Halo2-H2 


18.22 


2.21 


101.5 


690 


132 


2.2 


Halo2-HD 


18.74 


1.93 


94.4 


647 


61 


12.5 


Halo3-H2 


17.05 


2.29 


109.3 


662 


132 


2.1 


Halo3-HD 


15.37 


4.57 


151.8 


953 


58 


9.4 



Note. — See Table [2] for a description of the table format. 



